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The action for a class of three-dimensional dilaton-gravity theories with a negative cosmological 
constant can be recast in a Brans-Dicke type action, with its free u) parameter. These theories 
have static spherically symmetric black holes. Those with well formulated asymptotics are studied 
through a Hamiltonian formalism, and their thermodynamical properties are found out. The theories 
studied are general relativity (to —> oo), a dimensionally reduced cylindrical four-dimensional general 
relativity theory (to = 0), and a theory representing a class of theories (lu = —3). The Hamiltonian 
formalism is setup in three dimensions through foliations on the right region of the Carter-Penrose 
diagram, with the bifurcation 1-sphere as the left boundary, and anti-de Sitter infinity as the right 
boundary. The metric functions on the foliated hypersurfaces are the canonical coordinates. The 
Hamiltonian action is written, the Hamiltonian being a sum of constraints. One finds a new action 
which yields an unconstrained theory with one pair of canonical coordinates {M, Pm}, M being the 
mass parameter and Pm its conjugate momenta The resulting Hamiltonian is a sum of boundary 
terms only. A quantization of the theory is performed. The Schrodinger evolution operator is 
constructed, the trace is taken, and the partition function of the canonical ensemble is obtained. 
The black hole entropies differ, in general, from the usual quarter of the horizon area due to the 
dilaton. 

PACS numbers: 04.60.Ds, 04.20.Fy, 04.60. Gw, 04.60.Kz, 04.70.Dy 

I. INTRODUCTION 
A. Black hole thermodynamics 

Hawking radiation [3] is a phenomenon that emerges when one combines classical general relativity with quantum 
fields in a black hole background. It is thus a semiclassical phenomenon. The radiation is thermal and emitted at a 
given temperature T which depends on the black hole parameters. When this is linked to other black hole properties 
one finds that black holes indeed have entropy S and allow a thermodynamical description [2]. 



B. Path integral approach to black hole thermodynamics 



One can also get the thermodynamic properties of a black hole through a path integral approach to quantum 
gravity. In this approach one uses Feynman's idea that the amplitude < g2)^->2,t2\gi,Y,i,ti > to go to a bra vector 
state < g<ii S2 > *2 1 , with metric <?2 on a spatial hypersurface S2 at some generic prescribed time £2 , from one ket vector 
state, I (71 , Si, t\ >, with metric g\ on a spatial hypersurface Si at some generic prescribed time t±, can be found from 
the sum of all possible classical vacuum geometries, i.e., metric field configurations g, which take the prescribed values 
gi and gi on the surfaces Si and S2, respectively, i.e., < <?2j S2, ti\g\-, Si, t\ >= J D[g] exp (iS[g]), where D[g] is a 
measure on the space of all the gravitational field configurations, and S is the Lorentzian action of the theory. Now, 
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it is known that if one Euclideanizes the prescribed time t such that t = — z/3 and Euclideanizes the action / = — iS, 
one obtains a timeless amplitude, which is equivalent to the partition function Z = J D[g] exp (— S[g]). The partition 
function is then computed in the saddle point approximation, which finally connects to thermodynamics through the 
relation between Z and the Helmholtz free energy F, i.e., F = — TlnZ, where T = 1//3 is the temperature, and we 
put Boltzmann's constant /cb = 1- Since one has computed the partition function Z in this formalism one is working, 
a priori, in a well defined canonical ensemble, where the temperature T is given. However, solitary black holes are 
unstable against thermal fluctuations, indeed the heat capacity is in general negative. This means that the canonical 
ensemble is not defined at all and one should work instead in a microcanonical ensemble [jj. In the microcanonical 
ensemble the total energy E of the system, rather than the temperature T, is given, and the ensemble is usually 
well defined. On the other hand, calculations in the canonical ensemble are much easier than in the microcanonical. 
Thus, in order to stick to the canonical ensemble in a coherent fashion, one has to go beyond a solitary black hole, 
by enclosing it in an appropriate finite box. In this way, the path integral formalism and its inherent connection to 
the canonical ensemble, makes sense, if one carefully chooses the boundary conditions by fixing the local temperature 
at the surface of the box itself, as was done for the Schwarzschild case in [3] , and generalized to an arbitrary static 
field in []|. This can also be naturally extended to the grand canonical ensemble when one includes charge [f}]. If 
further, one puts a cosmological constant into the system, i.e., one deals with an asymptotic anti-de Sitter spacetime, 
one can dispense with the box, since the cosmological constant itself yields a natural box 0|. This approach, called 
the Brown- York formalism, has been further used a number of times in various different dimensions and for several 
different theories of gravity containing black holes. One should note that since the entropy and the other laws of black 
hole thermodynamics can be derived without considering matter fields, one finds that the black hole entropy is really 
an intrinsic entropy pertaining to the black hole geometry itself. 

C. Hamiltonian approach to black hole thermodynamics 

There are many other approaches to calculate the entropy and the thermodynamics of a black hole. The route we 
want to follow, motivated by the success of the path integral approach, is to build a Lorentzian Hamiltonian classical 
theory of the gravity in question, and then obtain a Lorentzian time evolution operator in the Schrodinger picture. 
Afterward one performs a Wick rotation from real to imaginary time, in order to find a well defined partition function. 
In more detail, the prescription implicit in this approach is: find the Hamiltonian of the system, calculate then the time 
evolution between a final state and an initial state, i.e., between the bra and ket vectors of those states, < #2, £2,*2| 
and |<?i,£i,ti >, and then Euclideanize time. Here, the amplitude to propagate to a configuration < (72, , ^2 1 from 
a configuration \g\, Si, ti > , is represented by < (72, £2, t%\ exp (—iHfa — ti)) \gi, £1, ii > in the Schrodinger picture. 
Euclideanizing time, t% — 1\ = —ifl and summing over a complete orthonormal basis of configurations g n one obtains 
the partition function Z = ^ exp (— (3E n ), of the field g at a temperature where E n is the eigenenergy of the 
eigenstate gn- This route is based on the Hamiltonian methods of [1, 0, [l(J, EH • It was developed by Louko and 
Whiting i n fl2ll fo r the specific problem of finding black hole entropies and thermodynamic properties, and further 
applied in |l3Lll4llT3.IT3| by Louko and collaborators, in [H by Bose and collaborators, and in [l§.[l9l| by Kunstatter 
and collaborators. This approach also points to an entropy where the degrees of freedom are in the gravitational field 
itself, since nowhere one mentions matter fields. 

The Louko- Whiting method [H] relies heavily on the Hamiltonian approach of [ll|, which in turn is an important 
ramification of the Arnowitt, Deser, and Misner approach, the ADM approach [§] (see also 0|), when applied to the 
full vacuum Schwarzschild black hole spacetime. This spacetime is better described by a spherically symmetric white 
hole plus black hole plus two asymptotically flat regions. These regions are well pictured in a Kruskal, or perhaps 
better, in a Carter-Penrose diagram. Indeed, Kuchaf [lH in studying within this formalism the Schwarzschild black 
hole found the true dynamical degree of freedom of the phase space of such a spacetime, by considering the spacelike 
foliations of the full manifold. This degree of freedom is represented by one pair of canonical variables. This pair 
is composed of the mass M of the solution and its conjugate momentum, which physically represents the difference 
between the Killing times at right and left spatial infinities. In [12], this method was adapted by considering a spacelike 
foliation to the right of the future event horizon of the solution, enabling one to find the corresponding reduced phase 
space. From this one can obtain a Hamiltonian H and thus the Lorentzian time evolution operator exp(-iHt) in the 
Schrodinger picture. 

This method has been applied for various theories of gravity and in several different dimensions. In [12J, the first 
paper of the series, a vacuum Schwarzschild black hole in four dimensions in general relativity was placed in a box, 
defined as a rigid timelike frontier. The values of the metric functions were fixed both at the horizon and at the box. 
Then, in [l3[, the same procedure was applied to a vacuum dilatonic black hole in two dimensions. Here there was 
also a rigid frontier, where the values of the fields, dilatonic and gravitational, were fixed. Afterward, in fl4j ]. the 
method was applied to the Reissner-Nordstrom anti-de Sitter black hole in four-dimensional Einstein-Maxwell theory 
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with a negative cosmological constant. This time there was no rigid frontier, so the far away asymptotic properties 
were defined at anti-de Sitter infinity. Next, in [TH], the same formalism is applied to the case of topological black hole 
solutions of the equations of motion for the four-dimensional Einstein-Maxwell-anti-de Sitter theory, leaving however 
part of the Hamiltonian method implicit. Finally, it was shown in that the entropy could also be calculated via this 
Hamiltonian method for five-dimensional spherically symmetric solutions of a one parameter family included within 
Lovelock gravity theory, where the action is comprised of the Ricci scalar term and the four-dimensional Euler density, 
i.e. the Gauss-Bonnet term, multiplied by an undetermined coefficient. There were some other developments. In [TtJ 
the Brown- York approach and the Louko- Whiting approach are compared. The main differences lie in the choice 
of boundary conditions, resulting in the fact that the Hamiltonian of Brown- York is the internal energy, whereas 
the Louko- Whiting Hamiltonian is the Helmholtz free energy. In [H, Gil the Louko- Whiting method was applied 
with some modifications to generic two-dimensional dilaton-gravity theories. Here we want to use the Louko- Whiting 
method to study the entropy and thermodynamic properties of three-dimensional black hole solutions in gravity 
theories with a dilaton, a negative cosmological constant and a free parameter u, which can be generically comprised 
in a three-dimensional Brans-Dicke theory with a cosmological constant. 



D. Three-dimensional dilaton-gravity black holes and their Hamiltonian approach to thermodynamics 

There is great interest in stud ying three-dimensional theories of gravity. One of the main reasons is richness of 
structure with one dimension less [2(1 I2H I22I I23I I24I I25L I26L I27L l28| . In fact, it is simpler to deal with certain concepts 
such as temperature, entropy, and flux of radiation when one works in three dimensions, and it also seems easier to 
try a three-dimensional quantum description of the black hole system. Theories in three dimensions, such as general 
relativity, have a rich structure. Indeed, three-dimensional general relativity with a negative cosmological constant has 
a black hole solution, the Banados-Teitelboim-Zanelli (BTZ) black hole (2Q.[2l|. The BTZ black hole is also a solution 
of string theory with a dilaton and other fields [22] , and there is also a three dimensional black string solution which can 
be recovered from an exact conformal field theory [23| . The existence of classical solutions depends on there being a 
negative cosmological constant. A negative constant is part of the three-dimensional black hole spacetimes structure, 
where the term in the action assumes a form given by —2 A 2 . The cosmological constant term is therefore always 
negative, and thus, the black hole solutions of three-dimensional theories are usually asymptotically anti-de Sitter. 
Beyond three-dimensional general relativity and effective three-dimensional string gravity theories there are other 
interesting three-dimensional theories with a cosmological constant. One of these can be recovered from dimensional 
reduction of four-dimensional cylindrical general relativity, yielding a three-dimensional gravity with a dilaton with a 
particular coupling for the kinetic dilaton term (24)]. It is then natural to set up a general three-dimensional dilaton- 
gravity theory, by including an uu paraneter, which yields different couplings for the kinetic term of the dilaton. This 
theory is then Brans-Dicke theory in three dimensions (25|. We will study in this connection the three-dimensional 
black hole solutions of [25], which are parametrized by the Brans-Dicke parameter u, where, in principle, u> can take 
any value, i.e., 00 > uj > —00. Since we want to perform a canonical Hamiltonian analysis, using an ADM formalism 
supplied with proper boundary conditions, it is necessary to pick up from the maze of solutions found in [25], only 
those that fulfill the boundary conditions we want to impose. First, we are interested only in solutions with horizons. 
Second, we want only solutions that are asymptotically anti-de Sitter at infinity. The cases of interest to be studied 
are then black holes for which uj — > ±00, 00 > uj > — 1, and — | > w > — 00. As in [25[ we choose three typical 
amenable cases where an analytical study can be done. These are u> — > 00 (or equivalently uj — > —00), uj = 0, and 
uj = —3. The theory for which u) — > 00 is general relativity, and the solution is the BTZ black hole [20(|. The theory 
for which uj = is equivalent to cylindrical four-dimensional general relativity and the corresponding black hole was 
found in The theory for which uj = — 3 is just a case of 3D Brans-Dicke theory, with a black hole solution that 
can be analyzed in this context [H]- If a quantum theory only makes sense if its classical form can be quantized by 
Hamiltonian methods, one should pick up only solutions which can be put consistently in a Hamiltonian form. 

Thus, using the prescription for the canonical variables, we foliate, following ADM, the three-dimensional spacetime 
with spacelike hypersurfaces of equal time, whose right end is at the anti-de Sitter infinity, and whose left end is at 
the bifurcation circle of a non degenerate Killing horizon. These canonical variables are then transformed to another 
set of appropriate canonical variables, directly related to the physical parameters of the black hole classical solutions. 
In turn, these latter variables are reduced to the true degrees of freedom, which can be represented by one pair of 
new canonical variables, the mass M and its conjugate momentum Pm- All this construction is limited to the right 
static region of the Carter-Penrose diagram of spacetime. This reduced theory is then quantized canonically. From 
the quantum theory obtained, one builds a time evolution operator in the Schrodinger picture. Upon continuation of 
this operator to imaginary time, inserting it into two generic base state vectors and taking the trace one obtains the 
canonical partition function, appropriate to a canonical ensemble. Note one needs to pay attention to the boundary 
conditions of the ensemble when trying to build an appropriate one for the geometries of a quantum theory of gravity. 
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These have to respect the stability properties of the semiclassical approximations of this quantum theory. Here 
these suitable boundary conditions are ensured by the fact that the black hole solutions are asymptotically anti-de 
Sitter. Now, the partition function obtained, given the right conditions, is dominated by classical Euclidean solutions. 
Following fl_4j we choose a renormalized Hawking temperature, which is taken as a fixed quantity in the canonical 
ensemble, due to the fact that the Hawking temperature of the black hole goes to zero at infinity. The entropy of the 
system can then be found using this whole formalism. We note that the entropy of black holes in three-dimensional 
dilaton-gravity theories of the kind we study, have been worked out in [26] through a completely different approach, 
namely, using known properties of two-dimensional conformal field theories. As well, the entropy and thermodynamics 
of the BTZ black hole have been exhaustively studied using a number of different methods, see [13] for a first study 
and [HI for a review. 

E. Structure of the paper 

The structure of the paper is as follows. In Sec. Ill Al we present the classical solutions of the three-dimensional 
dilatonic black holes, whose quantization through Hamiltonian methods we will perform. There is a free parameter 
w for which we choose three different values, u> — oo, 0, —3, corresponding to the BTZ black hole, the dimensionally 
reduce four-dimensional cylindrical black hole, and a three-dimensional dilatonic black hole, respectively. In Sec. Ill Bl 
we introduce the spacetime foliation through which we will define the canonical coordinates and which will allow us 
to write the action as a sum of constraints multiplied by their respective Lagrange multipliers. Then follow three 
sections, Sees. IIII[ HVj and El where we develop the thermodynamic Hamiltonian formalism for the to — oo, 0, —3 
black holes, respectively. In each section, we first give the metric and the dilaton fields (subsection A). Then we 
replace in the action of the theory the foliation ansatz, with the canonical coordinates, and determine the conjugate 
momenta of the coordinates. From there we determine the form of the Hamiltonian action through a Legendre 
transformation. It is then necessary to guarantee that there is a well defined variational principle, so as to allow the 
derivation of the equations of motion of the coordinates and respective momenta. This entails adding surface terms 
to the original action, which in turn forces the definition of asymptotic properties for the coordinates, momenta, 
and Lagrange multipliers. It is also necessary to fix some variables at the boundaries in order to have a well defined 
variational principle. After this is accomplished, we write the usual metric functions of a three-dimensional spherically 
symmetric, static black hole spacetime as functions of the canonical variables (subsection B). We then choose a new 
set of canonical variables and perform a canonical transformation. This new set has a physical meaning, i.e., one is the 
mass parameter of the black hole, and the other, the momentum conjugate to the mass, is the spatial derivative of the 
Killing time. We rewrite the action with the new variables and obtain a new pair of constraints with a redefined pair 
of Lagrange multipliers (subsection C). The new action is then reduced to its true degrees of freedom, with one pair of 
canonical variables, the mass and its conjugate momentum. The reduced Hamiltonian is nothing more than the initial 
surface terms written as functions of the new variables (subsection D). We quantize the theory canonically, defining 
a Hilbert space and an inner product therein, with the Hamiltonian now being an operator in this Hilbert space. 
We write the time evolution operator in the Schrodinger picture (subsection E). In order to obtain a thermodynamic 
description of the system we build a partition function in the canonical ensemble with fixed temperature by continuing 
the time evolution operator to imaginary time, and taking the trace in the mass eigenstates. This requires that we 
normalize the trace. Finally, through saddle point methods we obtain the thermodynamic functions (subsection F). 
In Sec. |Vl]we discuss and conclude. Throughout the paper we make h= 1, = 1, c= 1, and G = I. 

II. THE 3D BLACK HOLE SOLUTIONS WHICH ALLOW A PROPER THERMODYNAMIC 

HAMILTONIAN DESCRIPTION 

A. The 3D black hole solutions 

Three-dimensional black hole theories have been studied in, e.g., [13, E0, 113, 0, US EE H3, HI- A general 
action that incorporates most of these black holes is a Brans-Dicke action, with gravitational and dilaton fields and 
a cosmological constant. It is given by [25[ 

S=^- [ d^/^e' 2 * OR ~ Aoj{d(f) 2 + 4A 2 ) + B , (1) 
2tt J 

where g is the determinant of the three-dimensional metric g^, R is the curvature scalar, 4> is a scalar dilaton field, 
A is the cosmological constant, tu is the Brans-Dicke parameter, and B is a generic surface term. 
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The general solution for a static spherically symmetric metric, i.e., 
dimensions), is [25] 



circularly symmetric (since we deal with tree 



ds 1 
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dT 2 
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R 2 dip 2 , 
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w = -l, 
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(3) 
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where T,R are Schwarzschild coordinates, a is a constant related to the cosmological constant (see below), and b is a 
constant of integration (see below) , and the general solution for (j> is given by 



1 



For the constant a one has 



2(w + l) 
= constant , 

2|(w + l)A| 



ln(a R), 



w=-l. 



y/\(u + 2)(2u + 3)\ ' 

1, 
0, 

The constant b is related to the ADM mass of the solutions by 



«^-2,--,-l, 



a 
a 



OJ 
OJ 



- 2 --i 

-1 . 



(6) 
(7) 

(«) 

(9) 
(10) 

(11) 
(12) 

(13) 

Since we want to perform a canonical Hamiltonian analysis, using an ADM formalism supplied with proper boundary 
conditions, it is necessary to pick up only those solutions that fulfill the conditions we want to impose. First, we 
are interested only in solutions with horizons, so we take b to be positive. Second, apart from a measure zero of 
solutions, all solutions have a non-zero |A|. This does not mean straight away that the solutions are asymptotically 
anti-de Sitter. Some have one type or another of singularities at infinity, which do not allow an imposition of 
proper boundary conditions. So, from [HI with the corresponding Carter-Penrose diagrams, we discard the following 
solutions: oj = — 1 which is simply the Minkowski solution of a low-energy limit of string theory, — 1 > oj > — § since 
it gives weird conical singularities at Carter-Penrose infinity, and u = — | since all the Carter-Penrose boundary is 
singular. Thus, the cases of interest to be studied are black holes for which oj — > oo, oo > oj > — 1, and — I > oj > — oo. 
For b positive these solutions have ADM mass M positive, so well defined horizons. As in [25] we choose three typical 
amenable cases where an analytical study can be done. These are oj — > oo (or equivalently oj — > — oo), oj = 0, and 
oj = —3. The theory for which oj — > oo is general relativity, and the solution is the BTZ black hole [13, The 
theory for which oj = is equivalent to cylindrical four-dimensional general relativity and the corresponding black 
hole was found in [24|]. The theory for which oj = —3 is just a case of 3D Brans-Dicke theory, with a particular form 
for the kinetic dilaton term, which has a black hole solution that can be analyzed in this context [25| . 
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B. ADM form of the metric 



The ansatz for the metric and dilaton fields with which we start our canonical analysis is given by 

ds 2 = -N(t, rfdt 2 + A(i, r) 2 {dr + N r (t, r)dt) 2 + R{t, r) 2 dip 2 , 



(14) 



e~ 2 * = {aR{t,r))^ 



(15) 



6 



This is the ADM ansatz Q for the metric of spherically symmetric solutions of the three-dimensional Brans-Dicke 
theory. In this we follow the basic formalism developed by Kuchaf [ll[ . The canonical coordinates R, A are functions 
of t and r, R — R(t,r), A = A(i, r). Now, r = is generically on the horizon as analyzed in but for our 
purposes r — represents the horizon bifurcation point of the Carter-Penrose diagram [H| (see also In 
three spacetime dimensions the point represents a circle. The coordinate r tends to oo as the coordinates themselves 
tend to infinity, and t is another time coordinate. The remaining functions are the lapse N = N(t,r) and shift 
functions N r = N r (t, r) and will play the role of Lagrange multipliers of the Hamiltonian of the theory. The canonical 
coordinates R = R(t,r), A = A(t,r) and the lapse function N = N(t,r) are taken to be positive. The angular 
coordinate is left untouched, due to spherical (i.e, circular) symmetry. The dilaton is a simple function of the radial 
canonical coordinate, and it can be traded directly by it through equation (fl5|) . as will be done below. The ansatz 
(fT4l) - lfT5l) is written in order to perform the foliation of spacetime into spacelike hypersurfaces, and thus separates the 
spatial part of the spacetime from the temporal part. Indeed, the canonical analysis requires the explicit separation 
of the time coordinate from the other space coordinates, and so in all expressions time is treated separately from the 
other coordinates. It breaks explicit, but not implicit, covariance of the three-dimensional Brans-Dicke theory. Such 
a split is necessary in order to perform the Hamiltonian analysis. The metric coefficients of the induced metric on the 
hypersurfaces become the canonical variables, and the momenta are determined in the usual way, by replacing the 
time derivatives of the canonical variables, the velocities. Then, using the Hamiltonian one builds a time evolution 
operator to construct an appropriate thermodynamic ensemble for the geometries of a quantum theory of gravity. 
Assuming that a quantum theory only makes sense if its classical form can be quantized by Hamiltonian methods, one 
should pick up only solutions which can be put consistently in a Hamiltonian form. Thus, in the following we perform 
a Hamiltonian analysis to extract the entropy and other thermodynamic properties in the three three-dimensional 
Brans-Dicke black holes mentioned above, those for which w = oo,0, —3. 



III. HAMILTONIAN THERMODYNAMICS OF THE BTZ BLACK HOLE (U = OO) 

A. The metric 

For ui — > oo, the three-dimensional Brans-Dicke theory reduces to three-dimensional general relativity [25[. Then 
the general metric and dilaton solutions, given in Eqs. ^ and ©, reduce to the following 

ds 2 = - [(a Rf - M) dT 2 + (flj ^5 m + ^ ' (16) 
e" 20 = 1, (17) 

with M = b and a = \/2|A|. This is the BTZ black hole solution. Next, in Fig. [H we have the Carter-Penrose diagram 
of the BTZ black hole, where i?h denotes the black hole horizon radius, i?h = yM/ (2A 2 ), R — is the radius of the 
causal singularity, and R — oo is the spatial infinity, and we have discriminated the static and dynamical regions by 
roman numerals, namely the static right I and left F, and the dynamical future II and past IF regions. 



\ R=0 






/ / R =R h 


r \f 


I 







FIG. 1: The Carter-Penrose Diagram for the U) = oo case. 
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B. Canonical formalism 

The action |T]) with ui — > oo, and 2A 2 = Z~ 2 , where Z is the AdS length, becomes, excluding surface terms, 

S[A, R, A, P; N, N r ] = J dt J dr j— 2iV _1 AP + 2N- 1 N r R'k + 2N- 1 {N r )'KR + 2N- x N r A'R 

-2Af" 1 (Af r ) 2 A'i?' + 2NA~ 2 A'R' - 2N\- 1 R" + 4N\ 2 AR} , (18) 

where ' means derivative with respect to time t and ' is the derivative with respect to r, and where all the explicit 
functional dependences are omitted. There are two ways to arrive at the simple form of Eq. (fT8)l . One is to replace 
the ADM ansatz for the metric and the field (j>, given in (fl"4)l - ifl"5")l . into the action U). One obtains an expanded 
action, where one has polynomials on the canonical coordinates and the Lagrange multipliers, plus their derivatives, 
with respect both to time t and to space r. Afterward, one integrates in the angular coordinate tp and obtains a 
two-dimensional action. Of the terms composing it some are total derivatives which can be discarded. This may 
prove itself to be cumbersome. The other, simpler, way to arrive at the form of the action 1(18)1 . is to dimensionally 
reduce, through the Killing angle coordinate (p, the action before replacing the ansatz (fT4)l - lfT5)l . The new action is 
then a two-dimensional integral. In this way the total derivatives are seen straight away, allowing one to discard all the 
unnecessary surface terms. After the reduction and integration one can replace again the ansatz. Depending on the 
situation we use three different letters containing the same information but with slightly different numerical values. 
Thus, a, A, and I are related by a 2 = 2X 2 = l~ 2 , where I is the AdS length. From the action lfl8)l one determines the 
canonical momenta, conjugate to A and R, respectively, 

Pa = -2A- 1 {i? - R'N r } , (19) 

P R = -2JV -1 | A — (AN r )'^ . (20) 

By performing a Legendre transformation, we obtain 

U = N |-~PrPa- 2 A" 2 A' R' + 2A- 1 R" - 4A 2 Ai?| + N r {P R R' - P' A A) 

= NH + N r H r . (21) 

Here additional surface terms have been ignored, as for now we are interested in the bulk terms only. The action in 
Hamiltonian form is then 

S[A, R, P A , Pr; N, N r ] = J dt J dr |p A A + P R R - NH - N r H r } . (22) 
The equations of motion are 

A = ~NP R + (N r A)', (23) 

R = -^NP A - N r R' , (24) 

P R = 4X 2 NA~ {2N'A- 1 Y + (N r P R y , (25) 
P A = 4X 2 NR-2N'A- 2 R' + N r {P A )' . (26) 

In order to have a well defined variational principle, we need to eliminate the surface terms of the original bulk action, 
which render the original action itself ill defined when one seeks a correct determination of the equations of motion 
through variational methods. These surface terms are eliminated through judicious choice of extra surface terms 
which should be added to the action. The action l(22)l has the following extra surface terms, after variation 

Surface terms = (~N r P R SR + N r ASP A + 2NA- 2 R'SA - 2NA- 1 SR' + 2N'A~ 1 SR) |~ . (27) 

In order to evaluate this expression, we need to know the asymptotic conditions of each of the above functions 
individually, which are functions of (t, r). 
Starting with the limit r — » we assume 

A(t,r) = A + O(r 2 ), (28) 
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R(t,r 
PA(t,r 
P R (t,r 

N(t,r 
N r (t,r 



= Ro + R 2 r 2 + (9(r 4 ) 

= 0(r 3 ), 

= 0(r), 

= A^)r + 0(r 3 ), 

= 0(r 3 ). 



With these conditions, we have for the surface terms at r 

Surface terms L_ r 







(29) 
(30) 
(31) 
(32) 
(33) 



(34) 



Note that there are time dependences on the left hand side of the falloff conditions and that there are no such 
dependences on the right hand side, in the lower orders of the expansion in r. This apparent discrepancy stems from 
the fact that there is in fact no time dependence in the lower orders of the majority of the functions, but there may 
still exist such dependence for higher orders. Nevertheless, terms such as Rq are functions, independent of (t, r), thus 
constant, but undetermined. Their variation makes sense, as we may still vary between different values for these 
constant functions. 



For r — > oo we have 



A(t,r) 


= Ir- 1 + l 3 r 1 {t)r- 3 + O oc {r- 5 ) , 


R(t,r) 


= r + fpity- 1 +0°°(r- 3 ), 


P\(t,r) 


= 0°°{r- 2 ) 


P R (t,r) 


= 0°°(r- 4 ), 


N(t,r) 


= iJ(t,r)'A(t,r)- 1 (JV- + (t)+0 00 (r- B )) ) 


N r {t,r) 


= 0°°(r- 2 ). 



Here, as usual, I 2 = 2A 2 . These conditions imply for the surface terms in the limit r — > oo 

Surface terms^^ = 2S(M + {t))N + . 
where M+(t) = 2(rj(t) + 2p(t)). So, the surface term added to J22j is 

S ds [A, R; N] = J dt (2i?oA r iA,7 1 - A>+M 4 

What is left after varying this last surface term and adding it to the varied initial action (see Eq. (|22| ) is 



dt(2R 5(N 1 Ag 



6N+M 4 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 



(41) 



(42) 



(43) 



We choose to fix A^iAq 1 on the horizon (r = 0) and N + at infinity. These choices make the surface variation ((43] 
disappear. The term NiAq 1 is the integrand of 



n a (h)n a {t 2 ) = -cosh / dtN^A^ 1 ^) 



t2 



(44) 



which is the rate of the boost suffered by the future unit normal to the constant t hypersurfaces defined at the 
bifurcation circle, i.e., at r — » 0, due to the evolution of the constant t hypersurfaces. By fixing the integrand we are 
fixing the rate of the boost, which allows us to control the metric singularity when r —> [1 21 ] . 



C. Reconstruction, canonical transformation, and action 

In order to reconstruct the mass and the time from the canonical data, which amounts to making a canonical 
transformation, we have to rewrite the form of the solutions of Eqs. (fl6l) - lfT7j) . We follow Kuchaf [ll| for this 
reconstruction. We concentrate our analysis on the right static region of the Carter-Penrose diagram. In the static 
region, we define F as 



F(R(t,r)) = (aR(t,r)) 2 -b, 



(45) 
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and make the following substitutions 

T = T(t,r), R = R(t,r), (46) 

into the solution Q16|l . getting 

ds 2 = —(FT 2 — F _1 7? 2 ) dt 2 + 2(-FT'f + F- 1 R'R)dtdr + (-F(T') 2 + F^R 2 ) dr 2 + R 2 dtp 2 . (47) 

This introduces the ADM foliation directly into the solutions. Comparing it with the ADM metric lfl4|) . written in 
another form as 

ds 2 = -(N 2 - A 2 (N r ) 2 ) dt 2 + 2A 2 N r dtdr + A 2 dr 2 + R 2 dip 2 , (48) 
we can write a set of three equations 

A 2 = -F{T') 2 + F- 1 {R 1 ) 2 , (49) 

A 2 iV r = -FT'f + F^R'R, (50) 

N 2 -A 2 (N r ) 2 = F^-F^R 2 . (51) 

The first two equations, Eqs. (|49|) and Eq. I|50p . give 



-FTt + F-WR 
-F(T') 2 + F~ 1 (R 1 ) 2 ' 1 ; 

This one solution, together with Eq. (@9|), give 

N = - R ' T - T ' R . (53) 

^-F(T') 2 + F- 1 (R') 2 

One can show that N(t,r) is positive (see [lH). Next, putting Eqs. lf52j) - ([53"ll . into the definition of the conjugate 
momentum of the canonical coordinate A, given in Eq. lfT9|) . one finds the spatial derivative of T(t, r) as a function 
of the canonical coordinates, i.e., 

— T'= ^F-'APa . (54) 

Later we will see that —T' = Pm, as it will be conjugate to a new canonical coordinate M. Following this procedure 
to the end, we may then find the form of the new coordinate M(t, r), as a function of t and r. First, we need to know 
the form of F as a function of the canonical pair A , R. For that, we replace back into Eq. (|49|) the definition, in Eq. 
(EU, of T', giving 

R'V fP^ 2 



f = UJ-(tJ- < 55 > 

Equating this form of F with Eq. lj45"]) , we obtain 

M = 2X 2 R 2 - F , (56) 

where F is given in Eq. (|55|) and a 2 = 2A 2 , see Eq. ([8]). We thus have found the form of the new canonical coordinate, 
M. It is now a straightforward calculation to determine the Poisson bracket of this variable with Pm — —T' and see 
that they are conjugate, thus making Eq. (f54|l the conjugate momentum of M, i.e., 

Pm = l F ~ lAp A ■ (57) 

It is now necessary to find out the other new canonical variable which commutes with M and Pm and which 
guarantees, with its conjugate momentum, that the transformation from A, R, to M and the new variable is canonical. 
Immediately is it seen that R commutes with M and Pm- It is then a candidate. It remains to be seen whether 
Pr also commutes with M and Pm- As with R, it is straightforward to see that Pr does not commute with M and 
Pm, as these contain powers of R in their definitions, and {R(t, r), P R (t,r*)} — S(r — r*). So rename the canonical 
variable R as R — R. We have then to find a new conjugate momentum to R which also commutes with M and Pm, 
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making the transformation from {A, R; Pa, Pr } — > {M, R; Pm, P R } a canonical one. The way to proceed is to look 
at the constraint H r , which is called in this formalism the super- momentum. This is the constraint which generates 
spatial diffeomorphisms in all variables. Its form, in the initial canonical coordinates, is H r = -AP| + PrR'. In 
this formulation, A is a spatial density and R is a spatial scalar. As the new variables, M and R, are spatial scalars, 
the generator of spatial diffeomorphisms is written as H r = PmM' + P R R' , regardless of the particular form of the 
canonical coordinate transformation. It is thus equating these two expressions of the super-momentum H r , with M 
and Pm written as functions of A, R and their respective momenta, that gives us the equation for the new P R . This 
results in 



P R = Pr- 2A 2 P~ 1 AP a P + F^A^PaR" - F" 1 A" 2 P A A'i? / 



-F^A^PiR'. 



(58) 



We have now all the canonical variables of the new set determined. For completeness and future use, we write the 
inverse transformation for A and Pa, 



PmF) 



A = ((R') 2 P _1 
Pa = 2FP M ((R'fF~ 1 
In summary, the canonical transformations are the following, 



P M F) 



(59) 
(60) 



R 

M 



P 



R 



R, 

2X 2 R 2 - F , 

P R - 2A 2 F~ 1 AP A P + F^A^PaR" - F- x A- 2 PaA'R! 

-F^A^P^R', 
1 



Pm = ^F^AP 



(61) 



It remains to be seen that this set of transformations is in fact canonical. In order to prove that the set of equalities 
in expression l[61j) is canonical we start with the equality 



PaSA + P R SR - P M SM - P R 5R = [SRln 



2R! + AP A 



2R' - APa 
+ S [ AP A + R' In 



2R' - APa 



2R' + AP A 



(62) 



We now integrate expression (|62|) in r, in the interval from r = to r = oo. The first term on the right hand side of 
Eq. ((62]) vanishes due to the falloff conditions (see Eqs. pS ]) -([33 ]l and Eqs. l[35 jl -(|40" ll ). We then obtain the following 
expression 



/>oo />oo 

/ dr {P A SA + PrSR) - dr (P M 8M + P R SR) = Slu [A, R, P A ] , 
Jo Jo 



(63) 



where Suj [A, i?, Pa] is a well defined functional, which is also an exact form. This equality shows that the difference 
between the Liouville form of {R, A; Pr, Pa} and the Liouville form of {R, M; P R , Pm} is an exact form, which 
implies that the transformation of variables given by the set of equations ((6TJ) is canonical. 

Armed with the certainty of the canonicity of the new variables, we can write the asymptotic form of the canonical 
variables and of the metric function F(t, r). These are, for i — > 



F(t,r 
R(t,r 
M (t, r 

P R&r 
P M (t,r 



AR\A- 2 r 2 



0(r 4 ), 
Po + P 2 r 2 + 0(r 4 ) , 
2\ 2 R 2 + (4A 2 P P 2 - 4P|Aq : 
0(f), 
O(r). 



(64) 
(65) 
(66) 
(67) 
(68) 



For r — > oo, we have 



F(t,r) = 2\ 2 r 2 -2(r 1 (t) + 2p(t)) + 0°°(r- 2 ). 



(69) 
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R(t,r) = r + (2A 2 )- 1 ( oW^ _1 +0 00 (r- 3 ), (70) 

M{t,r) = M + (t) + 0°°(r- 2 ), (71) 

i^(t,r) = 0°°(r- 4 ), (72) 

P M (t,r) - 0°°(r- 5 ), (73) 



where Af + (i) = 2(r](t) + 2p(t)), as seen before in the surface terms (see Eq. (|4"2]l). 

We are now almost ready to write the action with the new canonical variables. It is now necessary to determine 
the new Lagrange multipliers. In order to write the new constraints with the new Lagrange multipliers, we can use 
the identity given by the space derivative of M, 

M' = -A" 1 [R'H + ^-P A H r ) . (74) 



2 

Solving for H and making use of the inverse transformations of A and Pa, in Eqs. lf59|) and |60| . we get 

H = - M ^ lR/ + m ^, (75) 

(f-w-fpiy 

H r = P M M' + P R R'. (76) 

Following Kuchaf jOj]], the new set of constraints, totally equivalent to the old set H(t, r) = and H r (t, r) = outside 
the horizon points, is M'(t, r) = and Pd (t, r) — 0. By continuity, this also applies on the horizon, where F(t, r) = 0. 
So we can say that the equivalence is valid everywhere. The new Hamiltonian, the total sum of the constraints, can 
now be written as 

NH + N r H r = N M M' + N R P R . (77) 

In order to determine the new Lagrange multipliers, one has to write the left hand side of the previous equation, Eq. 
(|77ll . and replace the constraints on that side by their expressions as functions of the new canonical coordinates, spelt 
out in Eqs. l|75"Tl - lf76"j) . After manipulation, one gets 

N M = r+N r P M , (78) 

(F-^Rr-FP^y 

N R = ^ -+N r R'. (79) 



(F- X (R') 2 - FP- 



M J 

Using the inverse transformations Eqs. l|5"9"l) -l|6"0" |) . and the identity R = R, we can write the new multipliers as 
functions of the old variables 

N M = -NF^R'A- 1 + ^N r F- 1 AP A , (80) 

N R = -^NP A + N r R' , (81) 

allowing us determine its asymptotic conditions from the original conditions given above. These transformations are 
non-singular for r > 0. As before, for r — > 0, 

N M (t,r) - -^(tJAo^+OCr 2 ), (82) 

N R (t,r) = 0(r 4 ), (83) 

and for r — > oo we have 

N M (t,r) = -N+(t) + 0°°(r- 4 ), (84) 

N R (t,r) = O 00 ^- 1 ). (85) 



The conditions (|82j) - (|85|) show that the transformations in Eqs. (|80j) - (|8T|) are satisfactory in the case of r — > oo, but 
not for r — * 0. This is due to fact that in order to fix the Lagrange multipliers for r — > oo, as we are free to do, we 
fix N+{t), which we already do when adding the surface term 



- / dtN+M+ (86) 
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to the action, in order to obtain the equations of motion in the bulk, without surface terms. However, at r = 0, we 
see that fixing the multiplier N M to values independent of the canonical variables is not the same as fixing -/ViAq" 1 to 
values independent of the canonical variables. We need to rewrite the multiplier N M for the asymptotic regime r — > 
without affecting its behavior for r — > oo. In order to proceed we have to make one assumption, which is that the 
expression given in asymptotic condition of M(t, r), as r — ► 0, for the term of order zero, Mq = 2A 2 i?o(i) 2 , defines Rq 
as a function of Mq, and Rq is the horizon radius function, Rq = R^(Mq). Also, we assume that Mq > 0. With these 
assumptions, we are working in the domain of the classical solutions. We can immediately obtain that the variation 
of Rq is given in relation to the variation of Mq as 

SR = (2X)- 2 Rq 1 5M , (87) 

where, as defined above, a 2 = I 2 = 2A 2 . This expression will be used when we derive the equations of motion from 
the new action. We now define the new multiplier N M as 



N M = -N M 



[l-g)+ g R (r 2 R 2 ) 



where g(r) = 1 + 0(r 2 ) for r — > and g{r) = 0°°(r 5 ) for r — > oo. This new multiplier, function of the old multiplier 
N M , has as its properties for r — > oo 

N M (t,r) = N + (t) + oo (r- 4 ) 1 (89) 

and as its properties for r — > 

N M (t > r) = N M (t)+O(r 2 ), (90) 

where Nq 1 is given by 

Nfi 1 = \ 2 N 1 RqR2 1 Aq . (91) 
When the constraint M' — holds, the last expression is 

A> M = A^- 1 . (92) 

With this new constraint N M , fixing A^A^ 1 at r = or fixing N M is equivalent, there being no problems with N^, 



which is left as determined in Eq. lj79j) . 

The new action is now written as the sum of «Se, the bulk action, and Sos, the surface action, 



S 



The new equations of motion are now 



M, R, P M , P R ; N M , Ar R ] = J dt dr (p M M + P R R - N R P R + N M [(1 - g) + gR (r 2 R 2 Q ) '] M') 

dt (2RqN^ 



~ N+M+) . 


(93) 


M = 0, 


(94) 


R = A^ R , 


(95) 


p M = (N M y, 


(96) 


PR = 0, 


(97) 


M' = 0, 


(98) 


^R = 0, 


(99) 



where we understood N M to be a function of the new constraint, defined through Eq. (|88j) . The resulting boundary 
terms of the variation of this new action, Eq. (|93|) . are, first, terms proportional to SM and 5~R on the initial and 

final hypersurfaces, and, second, J dt (^RqSNq 1 — M+5N+J . Here we have used the expression in Eq. f87j) . The 

action in Eq. (|93| yields the equations of motion, Eqs. ([94"l) - lf99"ll . provided that we fix the initial and final values 
of the new canonical variables and that we also fix the values of and of N+ . Thanks to the redefinition of the 
Lagrange multiplier, from N M to N M , the fixation of those quantities, and N+, has the same meaning it had 
before the canonical transformations and the redefinition of N M . This same meaning is guaranteed through the use of 
our gauge freedom to choose the multipliers, and at the same time not fixing the boundary variations independently 
of the choice of Lagrange multipliers, which in turn allow us to have a well defined variational principle for the action. 
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D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion 
94])- (J99j) allow us to write M as an independent function of space, r, 



M(t,r)=m(t). (100) 
The reduced action, with the constraints taken into account, is then 

m,p m ;A> M ,A> + ] = J dtp m m-h, (101) 

where 

/•OO 

Pm= / drP M , (102) 
Jo 

and the reduced Hamiltonian, h, is now written as 

h(m; t) = -2R h N^ + N+m , (103) 

with being the horizon radius. We also have that m > 0, according to the assumptions made in the previous 
subsection. Thanks to the functions Nq 1 (t) and N+(i) the Hamiltonian h is an explicitly time dependent function. 
The variational principle associated with the reduced action, Eq. (| 101|1 . will fix the value of m on the initial and final 
hypersurfaces, or in the spirit of the classical analytical mechanics, the Hamiltonian principle fixes the initial and final 
values of the canonical coordinate. The equations of motion are 

m = 0, (104) 
Pm = 2N ( f I (4X 2 R h )- 1 ~N + . (105) 

The equation of motion for m, Eq. (|104[) . is understood as saying that m is, on a classical solution, equal to the mass 
parameter M of the solution, Eq. (flB")) . In order to interpret the other equation of motion, Eq. (|105[) . we have to 
recall that from Eq. lj57j) one has Pm = —T', where T is the Killing time. This, together with the definition of p m , 
given in Eq. (|102|) . yields 

Pm = To - T+ , (106) 

where To is the value of the Killing time at the left end of the hypersurface of a certain t, and T + is the Killing time 
at spatial infinity, the right end of the same hypersurface of t. As the hypersurface evolves in the spacetime of the 
black hole solution, the right hand side of Eq. (|105|) is equal to Tq — T + , 



E. Quantum theory and partition function 



The next step is to quantize the reduced Hamiltonian theory, by building the time evolution operator quantum 
mechanically and then obtaining a partition function through the analytic continuation of the same operator |12l|-[16]. 
The variable m is regarded here as a configuration variable. This variable satisfies the inequality m > 0. The wave 
functions will be of the form ^(m), with the inner product given by 

= / ^mfe, (107) 

J A 

where A is the domain of integration defined by m > and /u(m) is a smooth and positive weight factor for the 
integration measure. It is assumed that /i is a slow varying function, otherwise arbitrary. We are thus working in the 
Hilbert space defined as J4? :— L 2 (A;^,dm). 

The Hamiltonian operator, written as h(t), acts through pointwise multiplication by the function h(m;i), which 
on a function of our working Hilbert space reads 



h{t)ip(m) = h(m; t)V>(m) 



(108) 
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This Hamiltonian operator is an unbounded essentially self-adjoint operator. The corresponding time evolution 
operator in the same Hilbert space, which is unitary due to the fact that the Hamiltonian operator is self-adjoint, is 



K(t 2 ;h) = exp 



12 



dt' h(t') 



This operator acts also by pointwise multiplication in the Hilbert space. We now define 

C*2 



T 



e 



dtN + (t) , 
o (*)■ 



dtN, M 



(109) 

(110) 
(111) 



Using (| 103[) . I|109p . IjllOp . and f| 1 1 If) . we write the function K, which is in fact the action of the operator in the Hilbert 

space, as 



K (m; T, 6) = exp [-imT + 2 i R h Q] 



(112) 



This expression indicates that K(t 2 ; t%) depends on t\ and t% only through the functions T and 6. Thus, the operator 
corresponding to the function K can now be written as K(T; 8). The composition law in time K(t3] t2)K{t2\ t\) = 
K{tz;ti) can be regarded as a sum of the parameters T and 9, inside the operator K(T; Q). These parameters are 
evolutions parameters defined by the boundary conditions, i.e., T is the Killing time elapsed at right spatial infinity 
and 6 is the boost parameter elapsed at the bifurcation circle. 



F. Thermodynamics 

We can now build the partition function for this system. The path to follow is to continue the operator to imaginary 
time and take the trace over a complete orthogonal basis. Our classical thermodynamic situation consists of a three- 
dimensional spherically symmetric black hole, asymptotically anti-de Sitter, in thermal equilibrium with a bath of 
Hawking radiation. Ignoring back reaction from the radiation, the geometry is described by the solutions in Eq. 
(|16p. Thus, we consider a thermodynamic ensemble in which the temperature, or more appropriately here, the inverse 
temperature (5 is fixed. This characterizes a canonical ensemble, and the partition function Z{(3) arises naturally in 
such an ensemble. To analytically continue the Lorentzian solution we put T = —i/3, and 6 — 2ni, this latter choice 
based on the regularity of the classical Euclidean solution. We arrive then at the following expression for the partition 
function 

Z{0) = Tr \k(-i/3, -2ni)] . (113) 

From Eq. (|112|) this is realized as 

/>oc 

Z(j3)= I /idm exp [-/3m + 47ri? h ] (m|m) . (114) 
Jo 

Since (m|m) is equal to 5(0), one has to regularize Eq. (|1 14|1 . Following the procedure developed in the Louko- Whiting 
approach [12], this means regularizing and normalizing the operator K beforehand. This leads to 



2ren(/3)=AA / fi dm exp [-/3m + 47ri? h ] , (115) 

J A 



where M is a normalization factor and A is the domain of integration. Provided the weight factor fi is slowly varying 
compared to the exponential in Eq. (|115p . and using the fact that the horizon radius is a function of m, the 
integral in Eq. (| 1 15|) is convergent. Thermodynamically, m is analogous to the energy of the system. Changing 
integration variables, from m to i?^, where 



m=2YR^ 1 (116) 



the integral Eq. I|115p becomes 



(f3)=Af JidRuexpi-Q, (117) 
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where A! is new the domain of integration after changing variables, and the function /^(i?^), a kind of an effective 
action (see 01), is written as 



h(R h ) := 2(3X'R^-AnR h 



(118) 



The new domain of integration, A' , is defined by the inequality i?^ > 0. The new weight factor Jl includes the Jacobian 
of the transformation, which amounts to dm/dR^. Since the weight factor Ji is slowly varying, we can estimate the 
integral of Z re n(/3) by the saddle point approximation. For that we have to calculate the critical points, i.e., we have 
to find the values of R^ for which the first derivative of /♦(-Rjj) with respect to R^ is zero. It happens for only one 
value of the domain 

K = j»- ^ 

In order to find out the nature of the extremum we evaluate the second derivative at the extremum. One finds 

8 2 h 



(120) 



As our domain starts at R^ = 0, we have that the extremum located at R^ = Ri is a minimum. Evaluating the 



action /* at R^ one obtains 



By Taylor expanding the action in the exponential of the integral, Eq. (|117p . 



I^R h )= h (R+)+ — 



1 d 2 h 



(R h ) 2 + 0((R h ) 3 



(121) 



(122) 



we can separate the terms in such a way that we obtain the following expression for the renormalized partition function 



ZreniP) = exp [-J.(i^)] M J JidR h exp ^A 2 ^ 



(123) 



The Taylor expansion was up to second order, and evaluated at the critical point R^, which makes the first order 
term of the expansion of /^(i?^) disappear. The term which can be put outside the integral is the zero order term, 
which is the value of /♦(i?^) at the extremum R^. The term left inside the exponential, 2/3A 2 i?^, is minus the second 
order term in the Taylor expansion, where all the higher orders have been ignored, as this is a good approximation, 
provided the weight factor is slowly varying. Finally, we may write the renormalized partition function as 



-ren 



/2tt 2 



where P is given by 



J JidR h exp(2f3X 2 Rl 



(124) 



(125) 



This P is a slowly varying prefactor and this approximation is better as we move to higher values of In 
the domain of integration the dominating contribution comes from the vicinity of i?^ = i?^. Leaving the explicit 
dependence of the partition function on the Ry. we write the logarithm of Z Ten (f3) as 



ln(Z ren (/3)) = InP - 2f3\ 2 (R+) 2 + AttR+ 



(126) 



By ignoring the prefactor's logarithm, which closer to is less relevant, we are able to determine the value of m 
at the critical point, where we find that it corresponds to the value of the mass of the classical solutions of the black 
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holes given in Eq. (flB")) . Thus, when the critical point dominates the partition function, we have that the mean energy 
(E) is given by 

(E) = -^lnZ ren (/3) w 2A 2 (i?+) 2 = m + , (127) 

where m + is obtained from Eq. (|116|) evaluated at R^. Using Eq. I|119p and Eq. (j 1 16[> . we obtain the temperature 
of the black hole T = (3^ 1 

^r) ' (128) 

where m + is the function in (| 1 16)1 evaluated at i?^. Inverting Eq. (jT28j) , and writing the inverse temperature as /?, 
we write the function m + (/3) 

m+ (/ ? ) = ^2" ( 129 ) 

We see that dm + /d(3 < 0, which through C = —(3 2 (d (E) /d(3) tells us that the system is thermodynamically stable. 
The entropy is given by 

S = (l - (In Z ren (/3)) w 47ri?+ . (130) 

This is the entropy of the BTZ black hole (see [H] and also [IE EE HI)- 



IV. HAMILTONIAN THERMODYNAMICS OF THE GENERAL RELATIVISTIC CYLINDRICAL 
DIMENSIONALLY REDUCED BLACK HOLE (LO = 0) 

A. The metric 



For uj = 0, the corresponding three-dimensional Brans-Dicke theory is obtained from the cylindrical dimensionally 
reduced black hole of four-dimensional general relativity 0, HB|- Then general metric in Eq. <[2j) and the <p field m 
Eq. ([6]), reduce to the following 



ds 2 



', ,i M 
(aR) 



2{aR) 



dT " + i m 2 _ M + R2 V , ( 131 ) 

\ aU ) 2(aR) 



e-'f = aR, (132) 

with M = 2b and a = y^||A|. Unlike the BTZ solution, this solution l|13ip - (|132p has a metric function whose mass 

term depends on R. This behavior is similar to the Schwarzschild black hole metric function [24] . In Fig. [2] we have 
plotted the Carter-Penrose diagram of the black hole solution for u = 0, where again R^ is the horizon radius, given 
by the larger positive real root of i?^ = j^t^, R — is the radius of the curvature singularity, and i? = oo is the 
spatial infinity. The roman numerals mean the same as in Fig. [TJ 

B. Canonical formalism 

The action with uo = becomes, excluding surface terms, 

S[A, R, A, R; N, N r ] = J dt J dra j X 2 ANR 2 - N^RRA - ^N^AR 2 + N- 1 RR{AN r )' 

+N- 1 N r RR'A - N- 1 N r ARR'(N r )' - N' 1 (N r ) 2 RR' A' + N^^ARR' 
_IiY-i(Ar)2 A (i?') 2 - (A-^'RR'N - ^A-^R'fN - A _1 i£R"ivj , (133) 
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where ' means derivative with respect to time t and ' is the derivative with respect to r, and where all the explicit 
functional dependences are omitted. Depending on the situation we use four different letters containing the same 
information but with slightly different numerical values. Thus, a, A, I and a are related by a 2 = |A 2 = l~ 2 , and 
a = 4a, where I is, as usual, the AdS length. From this action one determines the canonical momenta, conjugate to 
A and R respectively 



P A = -aN^R^R- R'Ar^ , 
P R = -olN^IrIA- {AN r )'} - 
By performing a Legendre transformation, we obtain 



U = N \ -oc^R-^PaPr 



1 



l KR~ 2 Pl + aA- 



RR 



A\R-N r R 



aA- 2 RR'A' + -A-^R') 2 - \ 2 ARf 



(134) 
(135) 



+N r {P R R' - P A A} = NH + N r H r . 
The action in Hamiltonian form is then 

S[A, R, P A , P R - N, N r ] = [ dt 



dr 



{PaA 



P r R - NH - N r H r 



The equations of motion are 

Na- X AR- 2 P A 
-No^PaR- 1 + N r R' 

Na' 1 

- x R- 2 Pl - (Na)'RR'A- 2 



A = -Nar 1 Rr 1 PR 



R 



P R = -Nu^PaPrR- 2 



N r A 



APlRr 3 



P A = — Na 



- ((^'A- 1 ^)' 
^Na(R') 2 A- 2 + 



- (Na)(A- 1 R r )' + 
Na\ 2 R 2 + N r P A 



2N a\ AR 



(136) 
(137) 

(138) 
(139) 

(N r P R )' (140) 
(141) 



In order to have a well defined variational principle, we need to eliminate the surface terms of the original bulk action, 
which render the original action itself ill defined for a correct determination of the equations of motion through 
variational methods. Through the choice of added surface terms one can achieve this elimination. The action (|13T[) 
has the following extra surface terms, after variation 

. — lDJTD' i AF'A-lDtD ATT r>_ JT D I AfACn. i atdd'a— 2i 



Surface terms = a (-NA^RSR' + N'A RSR - N r P R SR + N r ASP A + NRR'A- 2 8A)\™ 



(142) 



In order to evaluate this expression, we need to know the asymptotic conditions of each of the functions of (t, r). 
Starting with the limit r-iOwe assume 



A(t,r) 


= A (<) + 


0{r 2 ), 


R(t, r) 


= Ro{t) + 


R 2 (t)r 2 + 0(r 4 ) 


P\{t,r) 


= 0(r 3 ) 




P R (t,r) 


= 0(r), 




N(t,r) 


= Ni{t)r 


+ 0(r 3 ) 


N r (t,r) 


= N[(t)r 


+ 0(r 3 ). 



(143) 
(144) 
(145) 
(146) 
(147) 
(148) 
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With these conditions, we have for the surface terms at r = 0, 

Surface terms| r=0 = — aNiRoA^ 1 5Rq 



In the same way, for r — > oo, 



A(t,r 
R(t,r 
P A (t,r 
P R (t,r 
N(t,r 
N r (t,r 



= Ir- 1 + l 3 r j {t)r" i + 0°°(r~ 5 ), 

= r + l 2 p{t)r-' 2 + 0°°(r- 3 ), 

= 0°°(r~ 2 ) 

= o°°(^ 4 ), 

= 0°°(r- 2 ). 



These conditions imply for the surface terms in the limit r — ► oo, 

Surface termsl,,^^ = a5{M + )N + , 
where M + (t) = ot(i](t) + 3p(t)). So, the surface term to be added to (|137|1 is 

5 aE [A, R;N} = J dt QatiZgJVkV - A>+Af 4 

What is left after varying this last surface term and adding it to the varied initial action (see Eq. (|137p ) is 

dt (-aRlSiNiAo 1 ) - 5N+M+ 



(149) 



(150) 
(151) 
(152) 
(153) 
(154) 
(155) 



(156) 



(157) 



(158) 



We choose to fix ATiA 1 on the horizon and N+ at infinity, which makes the surface variation (|158[) disappear. 



C. Reconstruction, canonical transformation, and action 

In order to reconstruct the mass and the time from the canonical data, which amounts to making a canonical 
transformation, we have to rewrite the general form of the solutions in Eqs. (|131[) - f|132[l . We again follow Kuchaf [ll| 
for the reconstruction. We concentrate our analysis in the right static region of the Carter-Penrose diagram. In the 
right static region we define F as 

F(R(t, r)) = (aR(t, r)f - , (159) 

aR(t, r) 

and make the following substitutions 

T = T{t,r), R = R(t,r), (160) 

in the solutions above, Eqs. f|131[l - (|132[l . getting 

ds 2 = -{FT 2 - F^R^dt 2 + 2{-FT'f + F- 1 R'R)dtdr + (~F(T') 2 + F" 1 R 2 ) dr 2 + R 2 d(p 2 . (161) 

This introduces the ADM foliation directly into the solutions. Comparing it with the ADM metric lfl4|) . written in 
another form as 

ds 2 = -(N 2 - A 2 (N r ) 2 ) dt 2 + 2A 2 N r dtdr + A 2 dr 2 + R 2 dip 2 , (162) 
we can write a set of three equations 

A 2 = -F(T') 2 + F- 1 (R') 2 , (163) 

A 2 N r = —FT'f + F~ X R'R , (164) 

N 2 -A 2 (N r ) 2 = FT 2 -F- 1 R 2 . (165) 
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The first two equations, Eq. (|163p and Eq. (|164|) . give 



-FT'T + F~ 1 R'R 

N = -F(T'y+F-\R>V ■ (166) 



This one solution, together with Eq. (|163p . give 



R'T ~ T'R , x 

N = ■ . 167 

y/-F(T') 2 +F- l {R'f V ; 

One can show that N(t, r) is positive. Next, putting Eq. (|166|1 and Eq. I|167p into the definition of the conjugate 
momentum of the canonical coordinate A, given in Eq. (j 134[> . one finds the spatial derivative of T(t,r) as a function 
of the canonical coordinates, i.e., 

-T' = a~ 1 R~ 1 F~ 1 AP A . (168) 

Later we will see that —T' = Pm, as it will be conjugate to a new canonical coordinate M. Following this procedure 
to the end, we may then find the form of the new coordinate M(t, r), also as a function of t and r. First, we need to 
know the form of F as a function of the canonical pair A , R. For that, we replace back into Eq. I|163p the definition 
of X" , giving 

R'V fP\^~ 



F =( a - Ui> ■ (169) 



Equating this form of F with Eq. I|159p , we obtain 



1 ~^%2 



M = -aR \^—R z -Fj, (170) 

where F is given in Eq. (j 169[) . We thus have found the form of the new canonical coordinate, M. It is now a 
straightforward calculation to determine the Poisson bracket of this variable with Pm = —T' and see that they are 
conjugate, thus making Eq. I|168p the conjugate momentum of M, i.e., 

P M = a^R^F^APA . (171) 

It is now necessary to find out the other canonical variable which commutes with M and Pm and which guarantees, 
plus its conjugate momentum, that the transformation from A, R to M and the new variable is canonical. Immediately 
is it seen that R commutes with M and Pm- It is then a candidate. It remains to be seen whether Pr also commutes 
with M and Pm- As with R, it is straightforward to see that Pr does not commute with M and Pm, as these 
contain powers of R in their definitions, and {R(t, r), Pr(<, r*)} = S(r — r*). So rename the new canonical variable 
R as R — R. We have to find a new conjugate momentum to R which also commutes with M and Pm, making 
the transformation from {A, R; Pa, Pr} — » {M, R; Pm, Pr } a canonical one. The way to proceed is to look at 
the constraint H r , which is called in this formalism the super-momentum. This is the constraint which generates 
spatial diffeomorphisms in all variables. Its form, in the initial canonical coordinates, is H r = Pr R' — AP A . In 
this formulation, A is a spatial density and R is a spatial scalar. As all new variables, M and R, are spatial scalars, 
the generator of spatial diffeomorphisms is written as H r = PrR' + PmM' , regardless of the particular form of the 
canonical coordinate transformation. It is thus equating these two expressions of the super-momentum H r , with M 
and Pm written as functions of A, R and their respective momenta, that gives us the equation for the new Pr. This 
results in 

P R = P R - ^-F^APaR - -R^APa + F^Pa^'A- 1 - F^A^PaA'R' + (R'^F^A^PaR- 1 

-F^A^P'aR' . {Ill) 

We have now all the canonical variables of the new set determined. For completeness and future use, we write the 
inverse transformation for A and Pa, 

A = ((R') 2 F _1 -PmFY , (173) 
Pa = aRPP M ((R') 2 P" 1 -P^P)"^ . (174) 
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In summary, the canonical transformations are 
R = R, 
M -- 



P 



R 



aR 



Pm = 



a 

To 
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-F^A^P^R', 
a^R^F^APA 



R 



R 2 - F 
F^APaR 



1 



R^APa + F^Pa^'A- 1 - F- l A- 2 P A A'R' + {R'^F^A^PaR^ 1 



(175) 



In order to prove that the set of equalities in expression (|175|1 is canonical we start with the equality 

aRR' + AP A 



PaSA + P R SR - P M 8M - P R 5R = -aRSRhi 



aRR' - AP A 



1 



+ 6 ( APa + -aRR' In 



aRR' 


-AP A 


) 


aRR' 


+ AP a 





(176) 



We now integrate expression (|176p in r, in the interval from r = to r — oo. The first term on the right hand side 
of Eq. ([176]) vanishes due to the falloff conditions (see Eqs. lfT43l) - lfT48l) and Eqs. JUnj-jUlJ). We then obtain the 
following expression 



/•oo />oo 

/ dr (P A SA + P R SR) dr (P M SM + P R SR) = 5u [A, R, P A ] , 

Jo Jo 



(177) 



where 8w [A, P, Pa] is a well denned functional, which is also an exact form. This equality shows that the difference 
between the Liouville form of {R, A; Pr, Pa} and the Liouville form of {R, M; P R , Pm} is an exact form, which 
implies that the transformation of variables given by the set of equations (| 1 T5|) is canonical. 

Armed with the certainty of the canonicity of the new variables, we can write the asymptotic form of the canonical 
variables and of the metric function F(t, r). These are, for r — > 



F(t,r 
R(t,r 

M(t,r 

P M {t,r 
For r — > oo, we have 



= 4P 2 (t)A (i)-V+O(r 4 ), 
= R (t) + R 2 (t)r 2 + (9(r 4 ), 

= ^a 2 R (t) 3 + L aRo (t)R 2 (t) (3a 2 R {t) - 64R 2 {t)A Q {t)- 2 ) r 2 + 0{j 



0(r), 
0(r) . 



F(t,r) 

R(t,r) 
M(t,r) 

Pu{t,r) 



r 2 -2{ V {t) + 2p(t))r- 1 +0°°(r- 2 ). 



a 
16 

r + 16p(t)a~ 2 r- 2 + 0°°(r" 3 ) 

0°°(r- 4 ), 
0°°(r- 6 ), 



(178) 
(179) 

(180) 

(181) 
(182) 



(183) 

(184) 
(185) 
(186) 
(187) 



where M + {t) = a(rj(t) + 3p(t)), as seen before in the surface terms (see Eq. (|157p ). 

We are now almost ready to write the action with the new canonical variables. It is now necessary to determine 
the new Lagrange multipliers. In order to write the new constraints with the new Lagrange multipliers, we can use 
the identity given by the space derivative of M, 

M' = -A' 1 (R'H + a _1 iJ- x P A Fr) • (188) 
Solving for H and making use of the inverse transformations of A and Pa, in Eqs. (| 1 T3|) and (jT74j) , we get 



= M'P^R + FPmPr 

{F-i(Ry-FP M y 

H r = P Af M' + P R R'. 



(189) 
(190) 
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Following Kuchaf fll|, the new set of constraints, totally equivalent to the old set H(t, r) = and H r (t, r) = outside 
the horizon, is M'(t,r) — and Po(t,r) — 0. By continuity, this also applies on the horizon, where F(t,r) — 0. So 
we can say that the equivalence is valid everywhere. So, the new Hamiltonian, the total sum of the constraints, can 
now be written as 

NH + N r H R = N M M' + N R P R . (191) 

In order to determine the new Lagrange multipliers, one has to write the left hand side of the previous equation, Eq. 
(|19ip . and replace the constraints on that side by their expressions as functions of the new canonical coordinates, 
spelt out in Eqs. (| 189[) - (j lQOp . After manipulation, one gets 

N M = NF^R' J . +N rp Mi (192) 

N R = NFP M r+ N r R'. (193) 

(F-^Y-FPlY 

Using the inverse transformations Eqs. I|173p - (|174p . and the identity R = R, we can write the new multipliers as 
functions of the old variables 

N M = -NF^R'K- 1 + a- 1 N r F- 1 R- 1 KP K , (194) 
iV R = -a~ 1 NR- 1 P A + N r R' , (195) 

allowing us to determine its asymptotic conditions from the original conditions given above. These transformations 
are non-singular for r > 0. As before, for r — ► 0, 

N M (t,r) = ~N 1 {t)A {t)R i (t)- 1 + O(r > ), (196) 

N R {t,r) = -2Nr(t)R 2 (t)r 2 +0(r 4 ) = 0(r 2 ), (197) 

and for r^oowe have 

N M (t,r) = -N + (t)+0°°(r- 4 ), (198) 

N R (t,r) = 0°°(r- 2 ). (199) 

These conditions (|196j) - (|199[) show that the transformations in Eqs. (|194[) - (|195[1 are satisfactory in the case of r — > oo, 
but not for r — * 0. This is due to fact that in order to fix the Lagrange multipliers for r — * oo, as we are free to do, 
we fix N + (t), which we already do when adding the surface term 



- / dtN+M+ (200) 

to the action, in order to obtain the equations of motion in the bulk, without surface terms. However, at r = 0, we 
see that fixing the multiplier N M to values independent of the canonical variables is not the same as fixing ./ViA^" 1 
to values independent of the canonical variables. We need to rewrite the multiplier N for the asymptotic regime 
r — > without affecting its behavior for r — * oo. In order to proceed we have to make one assumption, which is 
that the expression given in asymptotic condition of M(t, r), as r — * 0, for the term of order zero, Mo = -^a 3 Ro(t) 3 , 
defines Rq as a function of Mo, and Rq is the horizon radius function, Rq = Ry.(Mo). Also, we assume that Mq > 0. 
With these assumptions, we are working in the domain of the classical solutions. We can immediately obtain that the 
variation of Rq is given in relation to the variation of Mo as 

32 

SR Q = —a~ 3 R- 2 5M . (201) 

This expression will be used when we derive the equations of motion from the new action. We now define the new 
multiplier N M as 



N M = -N M 



(202) 
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where g(r) = 1 + 0{r 2 ) for r — > and <?(r) = 0°°(r 5 ) for r — > oo. This new multiplier, function of the old multiplier 



N M , has as its properties for r — > oo 



and as its properties for r — > 



where iV^ is given by 



7V M (£,r) = JV+(t) + 0°°(r- 5 ) : 
iV M (t,r)=7V M (t) + O(r 2 ), 



When the constraint M' = holds, the last expression is 

iV M = A^Aq 1 



(203) 
(204) 

(205) 
(206) 



With this new constraint N M , fixing AiAq 1 at r = or fixing N M is equivalent, there being no problems with N R , 
which is left as determined in Eq. (j 1Q3[) . 

The new action is now written as the sum of Ss, the bulk action, and Sdn, the surface action, 



M, R,P M ,P R ;N M ,N R 



dt 



dr [ P M M + P R R - N R P Ti + V 



R 



(I "//I 1-2///?n( ^-Rl 



dt ( -aR£N<" - N+M+ 



The new equations of motion are now 



M 

R = 

Pm = 

P R = 

M' = 

P R = 



= o, 

= N- 



R 



(N M y. 

o, 

o, 

0. 



M'j 
(207) 



(208) 
(209) 
(210) 
(211) 
(212) 
(213) 



Here we understood N M to be a function of the new constraint, defined through Eq. (|2Q2j) . The resulting boundary 
terms of the variation of this new action, Eq. (|207[1 . are, first, terms proportional to 5M and SR on the initial and final 

hypersurfaces, and, second, the term J dt ^aR^SNQ 1 — M + 8N + \ Here we have used the expression in Eq. (|201|) . 
The action in Eq. (|20T|1 yields the equations of motion, Eqs. (|208|) - (|213[> . provided that we fix the initial and final 
values of the new canonical variables and that we also fix the values of Nq ' 1 and of N + . Thanks to the redefinition of 
the Lagrange multiplier, from N M to N M , the fixation of those quantities, N^ 1 and N+, has the same meaning it had 
before the canonical transformations and the redefinition of N M . This same meaning is guaranteed through the use of 
our gauge freedom to choose the multipliers, and at the same time not fixing the boundary variations independently 
of the choice of Lagrange multipliers, which in turn allow us to have a well defined variational principle for the action. 



D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion 
(|208jl - (|213|) allow us to write M as an independent function of the radial coordinate r, 



M(t, r) = m(t) . 

The reduced action, with the constraints and Eq. (|214|) taken into account, is 



5* 



m,P m ; Jv o". JV + = / dt p m m h, 



(214) 



(215) 
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where 

Pm= / drP M , (216) 



and the reduced Hamiltonian, h, is now written as 

h(m;t) = -^aR^N* 1 + iV+m, (217) 



with being the horizon radius. We also have that m > 0, according to the assumptions made in the previous 
subsection. Thanks to the functions 2Vq (t) and N+(i) the Hamiltonian h is an explicitly time dependent function. 
The variational principle associated with the reduced action, Eq. l|215p . will fix the value of m on the initial and final 
hypersurfaces, or in the spirit of the classical analytical mechanics, the Hamiltonian principle fixes the initial and final 
values of the canonical coordinate. The equations of motion are 

m = 0, (218) 
32 

Pm = - J a- 2 K h 1 N™-N + . (219) 

The equation of motion for m, Eq. (|218|1 . is understood as saying that m is, on a classical solution, equal to the mass 
parameter M of the solutions in Eq. (|131[) . In order to interpret the other equation of motion, Eq. (|218|1 . we have to 
recall that from Eq. (|17ip one has Pm = —T', where T is the Killing time. This, together with the definition of p m , 
given in Eq. (|2 16[) . yields 

Pm = T - T+ , (220) 

where To is the value of the Killing time at the left end of the hypersurface of a certain t, and T + is the Killing time 
at spatial infinity, the right end of the same hypersurface of t. As the hypersurface evolves in the spacetime of the 
black hole solutions, the right hand side of Eq. (|219p is equal to To — T+. 

E. Quantum theory and partition function 

The next step is to quantize the reduced Hamiltonian theory, by building the time evolution operator quantum 
mechanically and then obtaining a partition function through the analytic continuation of the same operator |12l|-[16]. 
The variable m is regarded here as a configuration variable. This variable satisfies the inequality m > 0. The wave 
functions will be of the form ip(m), with the inner product 

(ip,x) = [ Vdmjjx, (221) 

J A 

where A is the domain of integration defined by m > and /u(m) is a smooth and positive weight factor for the 
integration measure. It is assumed that fi is a slow varying function, otherwise arbitrary. We are thus working in the 
Hilbert space defined as Jf := L 2 (A; [idm) . 

Again, the Hamiltonian operator, written as h(t), acts through pointwise multiplication by the function h(m;i), 
which on a function of our working Hilbert space reads 

h(t)V>(m) = h(m; t)ip(m) . (222) 

This Hamiltonian operator is an unbounded essentially self-adjoint operator. The corresponding time evolution 
operator in the same Hilbert space, which is unitary due to the fact that the Hamiltonian operator is self-adjoint, is 



K{t 2 ;ti) = exp 



-i 



dt' h(t') 



(223) 



This operator acts also by pointwise multiplication in the Hilbert space. We define again 



T := / dtN+(t), (224) 
Jti 

9 := j 2 dtN^(t). (225) 
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Using l|217p , l|223p , <|224[) , and l|225p , we write the function which is in fact the action of the operator in the Hilbert 

space 




This expression indicates that K(t2',ti) depends on t\ and t% only through the functions T and 8. We can now write 
the operator corresponding to the function K as K(T; 8). The composition law in time K{t^; t^Rfo; ti) = K(ts; t\) 
can now be regarded as a sum of the parameters inside K(T; 8). These parameters are evolution parameters defined 
by the boundary conditions, i.e., T is the Killing time elapsed at right spatial infinity and 8 is the boost parameter 
elapsed at the bifurcation circle. 



F. Thermodynamics 

We can now build the partition function for this system. The path to follow is to continue the operator to 
imaginary time and take the trace over a complete orthogonal base. Our classical thermodynamic situation consists 
of a three-dimensional circularly symmetric black hole with a particular dilaton function, asymptotically anti-de 
Sitter, in thermal equilibrium with a bath of Hawking radiation. Ignoring back reaction from the radiation, the 
geometry is described by the solution in Eqs. f|131[l - (|132[) . Thus, we consider a thermodynamic ensemble in which the 
temperature, or more appropriately here, the inverse temperature (3 is fixed. This characterizes a canonical ensemble, 
and the partition function 2(f3) arises naturally in such an ensemble. To analytically continue the Lorentzian solution 
to imaginary time, we put T = — if) and 8 = — 27ri, based on the regularity of the classical Euclidean solution. We 
arrive then at the following expression for the partition function 



Z(fi) = Tr K(-i/3, -2ni) 



From Eq. <|226[1 this is realized as 



fi dm exp 



-/3m + naR^ 



mm 



(227) 



(228) 



Since (m|m) is equal to (5(0), one has to regularize (|228[) . Again, following the Louko- Whiting procedure [l3|-[lB|, we 
have to regularize and normalize the operator K beforehand. This leads to 



-2ren(/3) = N I A* dm exp 
I A 



-[3m + iraR^ 



(229) 



where J\f is a normalization factor and A is the domain of integration. Provided the weight factor /1 is slowly varying 
compared to the exponential in Eq. <|229[1 . and using the fact that the horizon radius is a function of m, the 
integral in Eq. <|229[) is convergent. Thermodynamically, m is analogous to the energy of the system. Changing 
integration variables, from m to Ru, where 



m 



a R^ . 



(230) 



the integral Eq. I|229p becomes 



-rem 



(/?) = Af I JldR h exp(-J*) 



(231) 



where A' is the new domain of integration after changing variables, and the function /^(i?^), a kind of an effective 
action (see is written as 



(232) 



The domain of integration, A', is defined by the inequality i?^ > 0. The new weight factor Jl includes the Jacobian 
of the change of variables, which amounts to dm/dR^. Since the weight factor is slowly varying, we can estimate 
the integral of Z ren (/3) by the saddle point approximation. For that we have to calculate the critical points, which 
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in our case amounts to finding the values of for which the first derivative of I^iR-^) with respect to R^ is null. It 
happens for two different values of the domain 



o. 



Rt = 2 6 tt (3/3 a 2 ) . 

In order to find the local extrema we evaluate the second derivative at these two points. One finds 

d 2 h 



dRl 



d 2 h 



HT 



*h 



-2ira . 



— 2na . 



(233) 
(234) 



(235) 
(236) 



Our domain starts at i?^ = 0, which is a local maximum. The global extremum, which is a minimum, is located at 
= i?^. Evaluating the action /* at R^ one obtains 



h(R+) = 2- 5 (3a 3 (R+) 3 -ira(R+) 2 
Substituting Eq. (|234j) into Eq. <|237|) gives 

I*(Rt) = -2 12 3- 3 ir 3 f3- 2 a- 3 . 



(237) 



(238) 



From Eq. (|238p one sees that I*(R^) < 0. By Taylor expanding the action in the exponential of the integral, Eq. 
((2311), i-e-, 



h(R h )=h(R+) + 



dh 



dRy 



1 (PI* 



Ri 



(R h ) 2 + 0((R h ) 3 )., 



(239) 



we can separate the terms in such a way that we obtain the following expression for the renormalized partition function 



-Zren(/3) = exp ~h(Rp Af I fidR^exp 



-TraR-r 



(240) 



The Taylor expansion is up to second order, and evaluated at the critical point R^, which makes the first order term 
of the expansion of I^(R] 1 ) disappear. The term which can be put outside the integral is the zero order term, which 
is the value of /^(i?^) at the extremum R^. The term left inside the exponential, — naR 2 ^, is minus the second order 
term in the Taylor expansion, where all the higher orders have been ignored, as this is a good approximation, provided 
the weight factor is slowly varying. Finally, we may write the renormalized partition function as 



2renW = P eM^~^P~^~ 3 ] 



where P is given by 



P = N I n dRfr exp 



(241) 



(242) 



This P is a slowly varying prefactor and this approximation is better as we move to higher values of |/*(i?^)|. In 
the domain of integration the dominating contribution comes from the vicinity of R^ = R^. Leaving the explicit 
dependence of the partition function on the R^ we write the logarithm of Z Ten {(3) as 



ln(Z re -n(/3)) = lnP + 7ra(i?+)" - 2-"f3a 3 {R+) 



(243) 



By ignoring the prefactor's logarithm, which closer to is less relevant, we are able to determine the value of m at 
the critical point, where we find that it corresponds to the value of the mass of the classical solution of the black hole 
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given in Eq. I|13ip , Thus, when the critical point dominates the partition function, we have that the mean energy 
(E) is given by 

(E) = -^lnZ ren (/3) « 2" 5 a 3 (i?+) 3 = m+ , (244) 

where m + is obtained from Eq. (|230[) evaluated at i?^ . By replacing the value of in Eq. (|234[) into Eq. l|230j) , we 
write the temperature T = 



3* |A| 3 m-< 
2§tt 3 



(245) 



where m + is the function in (|230|1 evaluated at R^. Inverting Eq. (|245[) we obtain the function m + (/3), where (3 1 
is the inverse temperature, 

m+(/3) = 2 13 a- 3 7r 3 (3/?)- 3 . (246) 

We see that dm + /d/3 < 0, which through the heat capacity C = —0 2 {d{E) /d/3) tells us that the system is thermo- 
dynamically stable. The entropy is given by 

S=(l- (lnZrenf/?)) « ™R+ 2 = 4™i?+ 2 . (247) 

This is the entropy of the three-dimensional black hole, which is a solution of the three-dimensional dilaton-gravity 
theory obtained through dimensionally reduced cylindrically general relativity (see also [Ifj])- Note that the entropy 
obtained in Eq. I|247p is not proportional to the area (i.e, circumference in the three-dimensional case) as the entropy 
of the BTZ black hole discussed in Section Hill This is certainly due to the presence of the coupling of the dilaton to 
the metric. 

V. HAMILTONIAN THERMODYNAMICS OF A REPRESENTATIVE DILATONIC BLACK HOLE 

(u = -3) 

A. The metric 

For lo = —3, the general metric in Eq. (J2j) and the <f> field in Eq. (J6J reduce to the following 



{aR) z - 2Va R M dT z + = — + R dtp 1 , (248) 

- (a R) 2 — 2 v aRM 

e" 20 = -^=, (249) 
ya R 

with 2M = b and a = ^p|A|. In Fig. [3] the Carter-Penrose diagram of the black hole solution for the case uj = —3 is 
shown, which is analogous to the case uj = 0, where again the singularity at R = is a curvature singularity. Since 
we are now familiar with the whole formalism, we will be briefer in this section omitting several of the details. 

B. Canonical formalism 

With uj = — 3 the solution implies that e _2< ^ = (aR)~i. The value of u> also implies that a = ^p|A| (see Eq. JSJ). 
The action then becomes, up to surface terms, 

S[A, R, A, R; N, N r ] = j dt dr ^i(a R)~^N\ 2 AR — (a R)~^ N~ 1 AR + ^(a R)~^ N~ 1 AR~ 1 R 2 

+ {aR)-iN' 1 R(N r A)' + (a R)~v N^^AR' - (aR)~iN~ 1 N r (N r ) 'AR' 
-{aR)-^N- 1 (N r ) 2 A'R' - ^(a R)~^ N' 1 N r AR~ X RR' + ^(a R)~^ N-\N r ) 2 AR' 1 (R') 2 

+ ^(aR)-?NA- 1 R- 1 {R') 2 - (a R)~t N^yR' - (a R)-? NA^R'X . (250) 
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Depending on the situation we use three different letters containing the same information but with slightly different 
numerical values. Thus, a, A, and I are related by a 2 = -yA 2 = §£~ 2 , where I is defined now as the AdS length. From 
the action above, we obtain the conjugate momenta 



Pa 



R - R N 1 



-N- l (aR)-Ti (A 
p R = -N~\aR)-b |a- (AN r Y - ^AR-\R- R'N r )^ 



(251) 
(252) 



By performing a Legendre transformation we obtain the Hamiltonian, which is a sum of constraints, i.e., 



H = N\-{aR)?P K P R + ^a*R-*kPl + {aR)-? 



1 

—t 

4 

N r {P R R' + A{P A )'} = NH + N r H r 



(A- 1 )^' + A^R" - -A-\R') 2 R- 1 



4(aR)-i\ 2 AR 



We can now write the action in Hamiltonian form, which reads 



S[A, R, P A , Pr] N, N r ] = J dt dr |f a A + P R R - NH - N r H r } 



(253) 



(254) 



with the constraints defined in Eq. (|253[) . From here we derive the equations of motion for the canonical variables 
and respective canonical momenta 



A 
R 

Pa 

Pr 



-N(aR)?P R + ^Na^R-^AP A + (N r A)' , 
-N(aR)ip A + N r R' , 



(255) 
(256) 



1 



Na^R-^-Pl - (N(aR)-t ) i?'A" 2 - jN(aR)~^(R') 2 R~ x A~ 2 + AN(aRy^\ 2 R + N r P A , (257) 



1 



NaiR-iP A P R + ^Na^R'^APl - Na^^(A^ 1 )' (R~^)' 
8 

jN(aR)-iA- 1 R- 1 R" - ( (NaT?)' A~ l R~?\ - ^Na'? A- x (R-i)" 



-2Na-?\ 2 AR-? + (N r P R )' . 



(258) 



For a correct variational principle to be applied, we have to find out what surface terms are left over from the variation 
performed, with the purpose of deriving the equations of motion in Eqs. (|255p ~ (|258p . These surface terms are, 

Surface terms = N(a R)~? R' A~ 2 SA - N(aR)~^A^ x 6R' + ((a R)~^' A~ X 5R + ha R)~* A^R^R'SR 

- N r P R SR + N r ASP A \™ . (259) 
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In order to know the form of l|259p for r — > 0, we assume 



A(t,r) 


= A + O(r 2 ), 


R(t, r) 


= R + R 2 r 2 + 0{r 4 ) , 


PA(t,r) 


= 0(r 3 ) 


P R (t,r) 


= 0(r), 


N(t,r) 


= N^ty + Oir 3 ) 


N r (t,r) 


= 0(r 3 ). 



The surface terms become 



Surface terms| r=0 = —Ni (a R) 2 A 1 8Ro ■ 
The asymptotic conditions for r — > oo are assumed as 



A(t,r) 


= — Ir^ 1 + 2*l2rj(t)r~* H 


-0°°(r- 3 ), 


R(t,r) 


= r + 2i/ip(t)r-5 + 0°°( r 




PA(t,r) 


= 0°°(r- 2 ) 




P R (t,r) 


= 0°°(r- 4 ), 




N(t,r) 


= i?(t,r)'A(t,r)- 1 (iV + (t) + 


0°°(r- 5 )), 


N r (t,r) 


= 0°°(r- 2 ), 





where I is the AdS length. The surface terms, Eq. I|259p . for r-*oo, are written as 

Surface terms^^^ = N+5M+ , 

where 

M+(i) = 2 5 3~^ 77(f) + 2 2 32 p(t) . 
Therefore, the surface term added to <|254[) is 

S ds [A,R;N] = fdt (2a- 1 (ai? )2Ar 1 A- 1 - JV+M+) 



(260) 
(261) 
(262) 
(263) 
(264) 
(265) 



(266) 



(267) 

(268) 
(269) 
(270) 
(271) 
(272) 



(273) 



(274) 



(275) 



With this surface term added we obtain a well defined variational principle. What remains after variation of the total 
action, Eq. ((254)) and Eq. 1(275)1 . is 



dt (2a-\a R )^ SiN^ 1 ) - SN+M+'j 



(276) 



The surface terms coming from the variation of the total action disappear as the result of the fixation of <5(A^A 1 ) 
on the horizon, r = 0, and of N + (t) at infinity, r — » oo. 



C. Reconstruction, canonical transformation, and action 

Repeating the steps of the two previous sections, we now give the main results concerning u> 
function is given by the expression 

F = (aRf - 2Va~RM . 
The Killing time T = T(t, r) is a function of (t, r), where we find that 

— T' = F- 1 AP A (aR)i , 



-3. The metric 



(277) 



(278) 
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which is equal to minus the conjugate momentum of the new variable M, or Pm — Pm(^ r) 
help of Eq. i|278p we find F = F(t, r) as a function of the canonical variables, 



F = {R') 2 A- 2 - aRPl 



-T'(t,r). With the 
(279) 



Summing up, the canonical transformations are 
R = R, 



M = haR)~* {(aR) 2 ~F) , 



P R - F' 1 R~ 1 A' 2 
1 



RR'kP' A 



(R') 2 P A A - RR"APa + RR'A'Pa 



-P^APa - 4A 2 P~ 1 P A PA . 



P M = F- 1 AP A (a J R)' 



(280) 



In addition, the relevant inverse transformations back to the old variable A and respective conjugate momentum Pa 
are 



^ -FPiy 



A = (F~ X {R 
P A = {aR)-?FP u {F- l {R'Y-FP\ 
We have to show that this transformation is canonical. This requires using the identity 

\(aR)-i& + AP A 



(281) 
(282) 



P A 5A + P R SR - P M SM - P R SR 



~(aR)-*6R\n 



1 



5 AP A + -{aR)-iR'\n 



(aR)sR' -AP A 

{aR)-^R' - AP A 



(aR)~iR' + AP A 



(283) 



We now integrate Eq. I|283p . in r, from r — to r = oo. The first term on the right hand side of Eq. (|283|1 vanishes 
due to the falloff conditions, Eqs. I|260p - (|265p and Eqs. (|26T|1 - (|2T2[i . We obtain then the following expression 



/>oo />oo 

/ dr (P A SA + PrSR) - dr (P M SM + P R SR) = 5u [A, R, P A ] , 
Jo Jo 



(284) 



where 5ui [A, R, Pa] is a well defined functional, which is also an exact form. As above, this equality shows that the 
difference between the Liouville form of {R, A; Pr, Pa} and the Liouville form of {R, M; P R , Pm} is an exact form, 
which implies that the set of transformations (|280|) is canonical. 

With this result, we write the asymptotic conditions of the new canonical coordinates for r — > 



F(t,r 
R(t,r 
M (i, r 

PM{t,r 



AR 2 2 A^ 2 r 2 + 0(r 4 



Ra + R 2 r 2 + 0(r 4 ), 
M + M 2 r 2 + 0{r 4 ). 
0(r), 
O(r), 



where we have 



2i-il- 2 Rl | ^P+Po 



2 



m 2 = -r 2 R R 2 1 ^r L R 



A-i, 

2 
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g 

— Z R0R2 — 4-R 2 ^o 

o 



(285) 
(286) 
(287) 
(288) 
(289) 



(290) 
(291) 
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For r — > oo, we have 

Fit, r) = ^r 2 r 2 - y f^ 1 ^ («»?(*) + 3*/>(*)) r* + 0~(r°) , (292) 

R(t,r) = r+ f^ r ^ +0 oo (r- 1 ), (293) 

M(t,r) = M + (i)+0°°(r-5), (294) 
P R (t,r) = 0°°(r- 4 ), (295) 

JV(t,r) = 0°°(r-*), (296) 

where M+(t) is defined in Eq. ((274|l . 

We now write the future constraint M' as a function of the older constraints 

M' = -A" 1 (r'H + (aR)?P A H r ^ . (297) 

Using the inverse transformations of A and Pa in Eqs. I|28ip and (|282|) . we obtain the same form for the old constraints 
as functions of the new variables 

H = (298) 

H r = P M M' + P R R'. (299) 
The new Hamiltonian, the total sum of the constraints, can now be written as 

NH + N r H R = N M M' + N R P R . (300) 
The new multipliers are, using Eqs. (|298p - (|300[) . 

N M = -+N r P M , (301) 

{F-^Rr-FP^y 

N R = NFP m T+ N r R'. (302) 

(F-^Y-FPl.Y 

Using the inverse transformations Eqs. i|28ip - (|282p . and the identity R — R, we can write the new multipliers as 
functions of the old variables 

N M = -NF^R'A- 1 + ArF _1 AP A (aP)5 , (303) 
N R = -NP A (aR)i + N r R' . (304) 



For r^Owe have, 



N M (t,r) = -^N^AoR.z 1 + 0(r 2 ) , (305) 
N R (t,r) = 0(r 4 ), (306) 



and for r — > oo we have 



N M (t,r) = -7V + (t)+0°°(r- 2 ), (307) 
N R {t,r) = 0°°(r~s). (308) 

Again, for r — ► 0, fixing N M (t,r), which means fixing A^ 1 (t)A P^" 1 , is not equivalent to fixing A^Ap 1 . It is thus 
necessary to rewrite N M for r — > 0. So, assuming Mo as a function of Po allows one to define the horizon radius 
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Rq = RiAMq). We are thus working in the domain where Mq > 0, the domain of the classical black hole solution. 
The variation of Rq is given in terms of the variation of M in the expression 



Remember that a 



¥i 



SR = -a _1 (a-R)"5<5Mo. 
A|. The new multiplier N M is obtained from the old N M as 



N 



M 



-N 



M 



(1 - g) + -gfR, 1 



(309) 



(310) 



where I 2 = 2A 2 , and g(r) = 1 + 0(r 2 ) for r — > and g(r) — 0°°(r 5 ) for r — > oo. This new multiplier, function of the 
old multiplier, N M , has as its properties for r — > oo 



and as its properties for r — > 



where Nq is given by 



A> M (<,r) = N+(t) + 0°°(r- 2 ) 
A> M (i,r)=A> M (t) + O(r 2 ), 



2X Z N 1 R R- L A . 
When the constraint M' = holds, the last expression is 

iV M = A^A- 1 . 



(311) 



(312) 



(313) 



(314) 



With this new constraint N M , fixing A^Aq 1 or fixing N M is equivalent, there being no problem with iV , which is 



left as determined in Eq. (|308|1 . 

The new action is then, summing both the bulk and the surface terms, 



S 



M, R,P M ,P R ;N M ,N R 



dt j dr yP M M + P R R- N R P 1 

dt (2a- 1 {aR )^N^ 1 - N+M+ 



R 



N 



M 



M' 



(315) 



The new equations of motion are now 



M = 


o, 


R = 




Pm = 


(N M ) 


P R = 


o, 


M' = 


o, 


P R - 


0. 



(316) 
(317) 
(318) 
(319) 
(320) 
(321) 

Here we understood N M to be a function of the new constraint, defined through Eq. I|310p . The resulting boundary 
terms of the variation of this new action, Eq. (|315|1 . are, first, terms proportional to SM and SR on the initial and final 

hypersurfaces, and, second, the term J dt ^2a _1 (a Rq^SNq 1 — M+SN+j. Here we used the expression in Eq. (|309p . 
The action in Eq. (|315p yields the equations of motion, Eqs. (|316|) - (|321[> . provided that we fix the initial and final 
values of the new canonical variables and that we also fix the values of Nq ' 1 and of N + . Thanks to the redefinition of 
the Lagrange multiplier, from N M to N M , the fixation of those quantities, Nq and N + , has the same meaning it had 
before the canonical transformations and the redefinition of N M . This same meaning is guaranteed through the use of 
our gauge freedom to choose the multipliers, and at the same time not fixing the boundary variations independently 
of the choice of Lagrange multipliers, which in turn allow us to have a well defined variational principle for the action. 
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D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion, 
Eqs. (|316p - <|321[) . allow us to write M as an independent function of the radial coordinate r, 



M (t, r) = m(i) . 

The reduced action, with the constraints and Eq (|322p taken into account, is 



m,P m ;^ M , w + = / <ftp m m-h, 



where 



Pm 



dr Pm ■ 



and the reduced Hamiltonian, h, is now written as 

h(m;t) = -2a' 1 (a R h )^N^ + N+m, 
with being the horizon radius. We also have that m > 0. The equations of motion are then 



0. 



Pm 



'(R^N™ - N + 



(322) 
(323) 

(324) 

(325) 

(326) 
(327) 



Here m is equal to the mass parameter M of the classical solution in Eq. (|248|) . The second equation of motion, 
Eq. I|327p . describes the time evolution of the difference of the Killing times on the horizon and at infinity, due to 
p m = T -T+ andEq. ll324l) . 



E. Quantum theory and partition function 

The steps developed in Sections IIIII and IIVI can be readily used here. So we do nont spell out the corresponding 
calculations in detail. 



F. Thermodynamics 



We can now build the partition function for this system, with the Hamiltonian given in Eq. (|325[) . The steps are 
the same as was the case with u — oo and w = 0. The thermodynamic ensemble is also the canonical ensemble. Thus, 
the operator is 



K (m;T, 9) = cxp imT + 2 i a^ 1 (a R h )^Q . 
Again, with T = — i/3 and 9 = — 2ni, we write the general form of the partition function as 



Z = Tr 



K(-i/3,-2ni) 



This is realized as 



Zren(P) =N f JldR h exp(-J*) 

J A' 



(328) 



(329) 



(330) 



where N is a normalization factor, A' is the domain of integration, and the function I*{R^), a kind of an effective 
action [4|. ^(R^), is written as 



I*(R h ) ■= f M h )i - ATra-^aR^ . 



(331) 
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The domain of integration, given by A', is defined by the inequality > 0. The new weight factor Jl includes the 
Jacobian of the change of variables, which amounts to dm/dR^. Its critical point is at R^ = |a _2 7r/3 _1 . The action 
evaluated at the critical point is 



J,(J2f) = -a-*/r* v 

From Eq. Q332p one sees that /^(i?^) < 0. It is seen that the critical point is a minimum. This implies that 

-2xen(/?) ~ -Pexp 
= Pexp 



(332) 



-§(a*J)* 



(333) 



We can now derive the basic thermodynamic results, as long as the approximation for the saddle point is valid, which 
means that we have to work in the neighborhood of the classical solution l|248p . where the critical point dominates. 
The expected value of the energy E is 



(E) 



d_ 

dp 



Here we have 



34 lApm" 1 



22 7T2 



m+(/3) 



(334) 



(335) 



(336) 



The derivative of m + (/3) with respect to (3 is negative, which means that the heat capacity is positive. The system is 
thus stable. Finally the entropy is given by 



S=(l- P^j (lnZ ren (/3)) « ^^^R+ . 
The result recovers the entropy for the three-dimensional dilatonic black hole with oj — — 3 (see also [If 



(337) 



VI. CONCLUSIONS 



We have continued the Louko- Whiting program of studying, through Hamiltonian methods, the thermodynamic 
properties of black holes in several theories in different dimensions (see and (l3l|-[l9|). 

Specifically we have calculated the thermodynamic properties of black hole solutions with asymptotic infinities 
that allow a well formulated Hamiltonian formalism in three-dimensional Brans-Dicke dilaton-gravity. Only certain 
values of the Brans-Dicke parameter ui are allowed in this juncture. The corresponding theories are general relativity, 
i.e., uj — > oo, a dimensionally reduced cylindrical four-dimensional general relativity theory, i.e., uj = 0, and a theory 
representing a class of dilaton-gravity theories, with a typical uj given by uj = —3. Within a three-dimensional context, 
we have built a framework for the classical Hamiltonian theory, where the metric functions are used as canonical 
variables, and where one foliates the spacetime with equal time spacelike hypersurfaces. These hypersurfaces go 
from the bifurcation circle of the horizon on the left, to the asymptotic anti-de Sitter infinity on the right. Then 
we have performed a canonical transformation based on a reconstruction of the canonical variables. This canonical 
transformation is an adaptation of Kuchaf's work to the Louko- Whiting method (l2l. HjjL fl4l. HR fl6| ■ In [ll| the 
boundaries are the left and right infinities of the Kruskal diagram (or if one wishes of the Carter-Penrose diagram), 
whereas in [H, [3, 0, [H, El and here, with a thermodynamic goal in mind, the boundaries are the bifurcation 
sphere, here a 1-sphere or circle, and the right infinity. When the classical equations hold one finds that the new 
canonical coordinate is indeed the physical parameter of the classical solutions of the dilatonic black holes in three 
dimensions, i.e., the black hole mass M. Its conjugate momentum, Pm = —T', is the spatial derivative of the Killing 
time. It survives the reduction, i.e., the elimination of constraints, becoming at last the difference between the 
Killing time at the bifurcation circle and at infinity. Its equation of motion is the equation for the evolution of the 
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difference between Killing times, with respect to the time parameter t of each hypersurface. The other coordinate, 
R, and respective momentum, Pj^, vanish after reduction, being pure gauge. With the new variables, come new 
constraints and new Lagrange multipliers. We then reduce the Hamiltonian to an unconstrained theory with one pair 
of canonical coordinates. Then one performs a quantization of the theory, where one replaces the functions of the 
physical observables for operators of a Hilbert space. Constructing then the Schrodinger evolution operator from the 
Hamiltonian and taking the trace on a suitable basis, we obtain the partition function of the canonical thermodynamic 
ensemble. This ensemble is well denned and under suitable conditions the classical Euclidean solutions dominate the 
partition function, yielding the thermodynamics of the systems. 

To sum up, as noted above, this formalism was previously applied to two [lj| (see also (l8l.fl9|). four [l2l [Ml. Uti [l7| . 
and five dimensions in several different theories. Here we have applied to three dimensions in a quite general 
dilaton-gravity theory. We have shown that in three-dimensional theories with well denned asymptotics the formalism 
fits well. As in other instances, the negative cosmological constant has had a stabilizing role to play here in what 
the thermodynamic results are concerned. Notwithstanding, several modifications were needed. First, in the powers 
of the fall-off conditions, and second due to the presence of the scalar dilaton field, which was reflected in the 
fact that it changed the powers of the radial R coordinate, to name a few. Although, in order to build a three- 
dimensional Lorentzian Hamiltonian theory, these modifications had to be made, we have derived a quantum theory 
and a statistical description of the systems in question, and found the corresponding thermodynamics, with precise 
values for the temperature and entropy of the black holes studied. 
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